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toward  the  calculation  of  the  influence  of  frontal  separations  for 

THE  SHORT  RANGE  PROGNOSES  OF  PRESSURE  AND  VERTICAL  VELOCITIES 


[This  is  a  translation  of  an  article  by  V.P#  Sa&okov,  In  Trudy,  TSIP, 

No  6j  1957 ,  pages  17-21 j  080:  4611-N] 

At  the  present  time ,  in  the  works  concerning  the  short-range 
prognosis  of  meteorological  elements  by  hydrodynamical  methods,  large-scale 
processes  are  studied.  The  usual  consideration  o£  such  processes  permits 
;;he  use  of  the  so  called  quasi-geostrophic  condition  for  the  analyses  and 
ution  of  the  hydrodynamics!  equations.  It  is  considered,  that  in  large- 
scale  process  the  chief  part  of  the  wind  velocity  is  determined  by  the 
geos trophic  relation. 

With  the  help  of  the  quasi-geos trophic  condition  it  is  possible  to 


write  down  the  basic  equations  of  dynamic  meteorology  (the  vorticity  equation 
and  adiabatic  equation)  approximately  in  the  form: 
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Here  the  pt/  p.  system  of  coordinates  is  taken: 

H  height  of  an  isobar ic  surface  /?*=-  constant, 
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^ y  —  the  Coriolis  parameter 

CO  angular  velocity  of  the  earth *s  rotation 

(jf)  geographical  latitude 


acceleration  of  gravity 

RX  T'  (  ifc  —  ,? ) 
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f-(—  gas  constant 
T  absolute  temperature 

^  y~~ ,  corresponding  adiabatic  and  ordinary  temperature  gradient 
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Vy1—  vertical  velocity  in  the  X ^  ccordiriate  system 
air  density 
velocity  components 
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The  system  of  equations  (1)~(2)  themselves  present  a  non-linear 

system  of  differential  equations  for  the  unknowns  /-/  and  *L^  .  Their 

solution  could  be  found  numerically  with  the  help  of  the  fast-operating 

numerical  calculating  machines.  In  general ^  system  (l)-(2)  has  an 

jaalytical  solution  for  the  functions  and  *c~*  under  certain 

Jtr 

.imitations* 

The  problem  concerning  the  calculation  of  frontal  separations  in  the 
rt  range  prognoses  is  set  down  by  the  following  method.  The  atmosphere  is 
considered  as  a  two  layered  fluid7  each  layer  of  which  corresponds  to  a  warn 
or  cold  air  mass*  In  each  layer  the  parameter  }T)  is  considered  as 
constant;  but  different  from  one  another. 


At  the  boundaries  let  us  use; 
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^  (the  earth  is  considered  as  flat), 
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where  Yet  is  the  normal  component  of  the  velocity  of  displacement  of 

the  front.  lad ex  1  corresponds  to  the  warm  air,  index  2  -  cold. 

Conditions  (5)  and  (6)  make  it  possible,  with  the  help  of  equation  (2), 
to  write  down  the  condition  for  the  function  *  at  the  surface  of 

separation.  Further,  equations.  (1)  and  (2)  could  be  solved  relative  to  *  %T 
Undar  the.  condition,  'that  the  parameter  Ol  **  constant.  Since  we  are 
considering  a  two-layer  atmosphere  with  a  constant  parameter  in  ea>.h  layer, 
therefore  it  is  possible  to  construct  an  equation  for  each  layer,  the 
analytical  solution  of  which  shall  serve  only  for  the  given  layer.  With 
the  help -of  the  conditions,  which  follow  from  (5)  and  (6),  both  solutions 
can  be  fitted  together. 

For  the  construction  of  the  solution  the  assumption  v?as  made,  whicn 
consists  In  that,  the  surface  of  separation  is  presented  as  a  quasi-horizontal 
surface.  As  a  result  we  arrive  at  the  necessity  of  solving  the  two-layer 
problem  with  the  supplementary  conditions  (5)  and  (6)  at  the  boundary  of  the 
two  mediums .  It  follows  that  we  keep  in  mind  for  this,  that  the  adopted 
assumption  cannot  be  interpreted  as  fact  by  the  exclusion  of  the  frontal  slope 
from  consideration,  nevertheless  we  shall  consider  completely  the  basic 
kinematie&l  and  dynamical  pecularities  of  the  front,  using  (5)  and  (6). 
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Here  th&  present  method  of  1  iasaticm.  1®  available*. 

Ih<s  solution  of  the  problea  is  sought  by  the  Fourier  aethod.  *icaUy 
th«  solutions  are  presented  la  ths  fossa  of  space  integrals:' 
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where  f  is  a  fcnc*?a  function*  which  depend*  on  the  characteristics  of 

the  thesaoSwric  field#  of  the  aenxMpfefere,  p  «■  -TT  —  height  of  the 

front,  J  - '  the  equation  of  the  frontal  surface 

in  she  eyfite®  of  coordinat**/ 
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are  function*  of  the  preblesB  considered*  ^  shall  not  give  the 

exprcfttiofift  for  fcfe*  Green8®  f motion ’ hare,  but  we  shall  only  note,  tfc&fc  they 
could  be  calculated  eaelytieally  or#  by  mams  of  the  numerical  solution  of 
the  differential  operator  of  the  ays tea  (!)“(*)• 


^Xfcs  solution  of  the  problem  in  such  a  form  was  first  obtained  by 
Euleeye?  sad  6.X.  M&x  chink  4 
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In  particular,  the  Greenes  functions  and  were 

computed  by  the  latter  method,  by  Kinkelmann  [2]. 

The  second  integrals  of  formulas  (7)  and  (8)  describe  the  chief 
effect,  arising  from  the  presence  of  the  frontal  division.  It  depends  on  the 
distribution  of  the  discontinuities  in  the  wind  field  at  the  frontal  surface 
and  from  the  slope  of  the  front.  Here,  integration  is  carried  out  along  the 
frontal  surface. 

Thus  we  considered,  in  a  complex  way,  a  model  of  the  atmosphere, 
which  includes  in  its  detail,  the  frontal  processes  with  the  usual  large  scale 
processes.  It  is  completely  obvious,  that  the  frontal  processes  appear  to 
attribute  to  the  processes  a  smaller  order  than  the  large  scale  effects.  This 
assertion could  be  confirmed  by  empirical  evidence. 

On  the  strength  of  the  above-stated,  it  is  possible  to  consider 
separately  the  large  scale  processes,  described  by  the  first  integrals  in  (7) 
and  (8),  and  the  frontal  processes  described  by  the  second  integrals. 

Independently  analyzing  the  contributions  of  the  two  integrals  on  the 
general  change  of  pressure,  it  follows  that  we  keep  in  mind,  that  the 
relation  between  the  characteristic  horizontal  and  time  scale  calls  for  a 
calculation  of  the  second  integral  for  smaller  time  intervals,  than  is 
usually  done. 

We  calculated  the  correlation  coefficient  between  the  actual  3  hour 
pressure  change  and  the  values  of  the  two  Integrals  calculated  for  that  same 
interval  of  time.  The  correlation  coefficient  was  found  equal  to  0.6  according 
to  the  results  of  9  cases  (over  100  points).  Hence  it  is  possible  to  conclude, 
that  the  frontal  terms  in  the  formulas • (7,8)  quite  satisfactorily  describes 
the  real  process  in  this  time  interval.  The  effect  of  the  large-scale 
process  in  this  case  appears  as  a  background,  on  which  the  effect  of  locals 


process  (in  particular  frontal  processes)  is  superimposed. 

Let  us  turn  to  consideration  of  the  vertical  motions.  We  shall  be 
interested  in  those  vertical  motions ^  which  arise  in  the  region  of  the  £ron<~ 
and'  connected  on  the  whole  with  the  discontinuities  in  the  wind  field.  Let 
us  note*  that  contemporary  methods  of  calculation  of  the  vertical  motions 
[l3>  1 3]  have  in  mind  those  motions ^  which  arise  as  a  result  of  large  scale 

exchange  between  the  kinetic^  potential  and  intenal  energies  of  the  air  mass 

in  the  atmosphere.  It  is  possible  to  conisdar  the- vertical  motions  of  such 
a  generation  as  a  background.  They  can  not  describe  those  details  in  the 
distribution  of  vertical  motions^  which  are  observed^  for  example^  ata  front 

i 

by  virtue  of  the  peculiarities  of  its  structure.  In  the  future,  we  shall 
only  consider  those  vertical  motions,  which  are  caused  by  the  presence  of  the 
frontal  separation,  not  dwelling  on  consideration  of  the  large  scale  vertical 
motions. 

Let  us  determine  the  vertical  velocity  from  equation  (1),  inta gratia 
It  with  respect  to  the  vertical  and  requiring  that  the  quantity  C 

becomes  zero  at  the  upper  and  lower  boundary  of  the  atmosphere.  Hot  paying 

attention  to  the  large  scale  motions,  described  in  equation  (1)  by  the 
terms,  contained  on  the  right  side  in  expressions  (7,8)  by  the  first 
integrals,  it  is  possible  to  write 
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I a  exprevsiaa*  (11,  12),  only  the  Greeks  functions  alone  depend  on 


expr 
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for  tsx&t  reason  the  integration  with 
and  differentiation  ( A  *'jyl  )  ■*«  be  carried 


out  only  ever  tba  Green1#  functions.  As  a  result  formulas  (9,10),  after 
substitution  of  egressions  (11,  tt>  into  them,  can  be  rewritten  in  the  form: 
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Ssre  Q-j.'’  C?J;>  are  ne»  Green*  s  functions.  These  functions 
could  be  tabulated  under  some  assumptions  [  %jt  ^  D  , 

Th&  (  Ue  ~~  tit.  /  ix  +!rC/‘ k  v.,_  j  w 

c&&  b$  supposing ^  th&t  £l  &  JL  *  &  2  ?C 

Trcrn  the  cKmditloa  &t  *&«  front 
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in  particular,  it  follows,  that 
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Tmz&f  va  m^ge  to  express  tfea  discontinuity  in  the  wind  at  the 
fto'it  mad  its  alopfi  through  the  p&eul&rifcies  of  the  baric  field  at  it. 

Xa  expressions  (13,14),  the  Grs*mfs  functions  possess  that  characteristic, 
that  they  das?.p  v&ry  rapidly  from ■  the  forecast  point,  this  wans  that  the 
^  surroundings  of  the  poiut  •  of  calculation  of  the  field  of  the  functions 
(  H$t  /i  )  b&va  a  weak  influence  on  the  rsnult  of  the.  calculation.  For 
that  rosaon  it  is  possible  to  simplify  formulas  <13,14),  combining  them  into 
one  by  th©  following  fans 

(IS)  *  ( f  t  %  a'  ( p*  JT ^  (  ^  ^  ^ 

XI  we  select  &r  tfce  unit  of  length  in  sp&cs  as  500  km  and  we  measure  f~f  in 
dkm  (decuas&ters),  to  d&t^rmi ne  all  the  meteorological  parameters  in  the  KTC 
.  syetess.  (meter^tea-eecoad),  then  for  the  value  y£?«  0.8  (2 loo.)  the  values  of 
the  coefficients  &  for  the  calculation  of  c"  in  ejo/12  hours  will  be 
equal  to:  £.  (0.8;  1.0)  «  0.16;  *t<0.8;  0.85)  *  4.68; 

(J.  (0.8;  0.7)  ■»  0.42;  #<0.8;  0.5)  »  0.28;  #  <0.8;  0,3)  *  0.06. 

The  calculation  of  vertical  ssotions  in  the  frontal  zone  roughly  comas  to  the 
followings 
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1.  The  location  of  the  front  is  found  on  all  the  /7er^ 

and  sAT^op  maps* 

2*  The  positions  of  the  front  at  all  levels  are  placed  on  a  separate- 
clean  form  of  a  synoptic  map  [as  if  we  project  the  fronts  at  all  standard 
levels  on  to  one  plane  (map)]. 

3*  The  necessary  number  of  points  are  marked  on  the  line  of  the  fron 
at  each  isqbaric  surface* 


4.  At  the  marked  points  the  Jacobian  /  />4/ 

_  2  Hi  w  ti  ///  ^  „  '2j££  * 

■  ^  J5?  r  ^  ^ 

•3  calculated.  The  calculations  of  the  derivatives  la  carried  out  in  finite 
differences  *  For  this  it  is  convenient  to  direct  the  //  axis  along  the 
normal  to  the  front  on  the  side  of  the  cold  air^  while  the  axis  is  alone. 


the  tangent  (figure  1). 
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we  obtain  . 
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5, .  All  the  values  (  )  9  calculated  at  the  points  of  the 

surface  front  were  then  multiplied  by  the  coefficient  &,  (1)  *  0.16„  The 

values  (Hi.}  Hi  J  corresponding  to  the  points  of  the  front  at  the  850  mb 

surface,  were  multiplied  by  CL  (0,85)  *^4,68  etc*  These  results  were 

written  down  at  the  corresponding  points  on  the  clean  blank  (form)^ 

is 

mentioned  in  (2)^  and  then  the  map  ± si  drawn. 

In  operational  work  it  is  possible  to  use  an  accelerated  method  of 
Iculation^  which  consists  in  the  determination  of  the  maximum  values  along 
cne  whole  front.  The  actual  or  prognostic  map  is  used  for  this 

purpose.  The  position  of  the  front  is  found  on  it  and  (  Hf  /  )  is 

calculated  at  the  points  of  interest*  Then  the  result  is  multiplied  by  the 
coefficient  equal  to  4. 


In  conclusion^  we  give  the  results  of  some  examples  of  the 

.jj( 

calculation  of  (■>  on  warm  fronts .  Calculations  at  cold  fronts  proved 

to  be  less  successful,  especially  in  the  -warm  time  of  the  year.  The 
calculation  of  C  in  mb/ 12  hr  presented  in  figures  2,  3,  4,  corresponds 
to  the  morning  of  March  24,  1953,  the  morning  of  September  14,  1S53  and  the 
evening  of  October  14,  1954  (  the  shaded  regions  are  clouds  at  the  time  of 
the  calcuhtions) . 
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